Cold atom gravimetry with a Bose-Einstein condensate 
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We present a cold atom gravimeter operating with a sample of Bose-condensed 87 Rb atoms. Using 
a Mach-Zehnder configuration with the two arms separated by a two-photon Bragg transition, we 
observe interference fringes with a visibility of (83 =b 6)% at T = 3 ms. We exploit large momentum 
transfer (LMT) beam splitting to increase the enclosed space-time area of the interferometer using 
higher-order Bragg transitions and Bloch oscillations. We also compare fringes from condensed and 
thermal sources, and observe a reduced visibility of (58 ±4)% for the thermal source. We suspect 
the loss in visibility is caused partly by wavefront aberrations, to which the thermal source is more 
susceptible due to its larger transverse momentum spread. Finally, we discuss briefly the potential 
advantages of using a coherent atomic source for LMT, and present a simple mean- field model to 
demonstrate that with currently available experimental parameters, interaction-induced dephasing 
will not limit the sensitivity of inertial measurements using freely- falling, coherent atomic sources. 



I. INTRODUCTION 



In principle, a light bulb and laser that have the same 
photon flux will yield the same precision in many shot 
noise limited optical measurements. In practice, it is of- 
ten the classical properties of optical lasers - the bright- 
ness, coherence, and low phase and amplitude noise - 
that enable a shot noise limited optical measurement at 
high flux. To date, inertial measurements using atom in- 
terferometers have primarily utilised cold thermal sources 
[T]-[3]. It is therefore of interest to investigate whether 
coherent, high brightness atomic sources such as Bose- 
Einstein condensates (BECs) will have similar advan- 
tages in such measurements. To date, the use of BEC 
has been largely neglected by the precision measurement 
community. One concern is that the comparatively high 
atom density will lead to interaction-induced dephasing, 
thereby limiting precision [4j [5] . 

Here we present results on a gravimeter based on 
atomic interference [6j [7] of Bose-condensed 87 Rb, and 
compare its performance to that achieved with a cold 
thermal sample in the same system. We observe an in- 
crease in fringe visibility when using a condensed source 
instead of a thermal one in an equivalent setup. With 
one of the most promising avenues for increasing sensitiv- 
ity being large momentum transfer (LMT) beamsplitting 
[8UTT] , we use LMT to increase our sensitivity to gravity, 
and maintain good fringe visibility for both Bragg- and 
Bloch-based LMT. Furthermore, using a simple model we 
demonstrate that dephasing in an expanded BEC will not 
limit the precision of inertial measurements. 
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II. BACKGROUND AND METHODS 

The operating principle of a gravimeter based on 
atomic interference has been described elsewhere [TJ [6] . 
Our gravimeter uses nth-order Bragg transitions [9j [12] 
H3] as our atom-optic beamsplitters (tt/2 pulses) and mir- 
rors (tt pulses) in a Mach-Zehnder (tt/2— tt— tt/2) configu- 
ration. These couple vertical momentum states separated 
by 2nHk, where k = |fc| is the wavenumber of the light 
and n an integer. For uniform acceleration the atomic 
phase evolution of each arm is identical, and the only in- 
terferometric phase contribution is from the atom-light 
interaction pQ: 

$ = n(fa - 20 2 + 03 ) = 2n k ■ gT 2 (1) 

where fa is the optical phase of the zth Bragg pulse, 
and T is the time between pulses. Scanning <E> results 
in fringes P — \{A + V cos<I>) in the relative population 
in state |po + 2nhk) : where we define V as the visibil- 
ity [33] and A the fringe offset. One infers <I> by oper- 
ating at mid-fringe, where the change in P for a given 
phase shift is maximal. For small shifts A<l> one obtains 
a signal of AP = = Vnk • AgT 2 . High signal 

gain thus requires having a high visibility and a large 
space-time enclosed area (ex 2nkT 2 ). BEC interferom- 
eters have already been shown to exhibit fringe visibil- 
ity close to 100% [15] . The evolution time T is typi- 
cally limited due to practical considerations such as the 
fall distance of the atoms [6;, and low- frequency phase- 
noise sensitivity [16]. Increasing the space-time area us- 
ing large-momentum-transfer (LMT) beamsplitting can 
be achieved using higher-order (n > 2) Bragg transitions 
[9] or Bloch oscillations [8| flO| fTTj . 

LMT techniques require the momentum width of the 
source along k to satisfy Ap <C hk : so that the entire 
cloud may be coupled to a single momentum state. As 
suggested in [17J, this requirement on Ap becomes more 
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stringent for increasing n in order to maintain a high 
LMT efficiency. This can be understood by considering 
the dispersion relation, AE n /Ap = p n /m = 2nhk/m, 
where E n is the energy of the nth Bragg resonance. This 
shows that a cloud of given momentum width Ap has a 
greater spread in energy at higher momentum. Thus, for 
a Bragg transition resonant with the centre of the cloud, 
there will be greater spread in detuning from resonance 
across the cloud for higher order transitions. 

Using Bragg spectroscopy [18] , we measure the l-cr mo- 
mentum width of our cloud to be Ap = O.lAhk after 
12ms of ballistic expansion (see Fig. [2^). In contrast, 
even the coldest reported thermal source at a tempera- 
ture of 150 nK has Ap = 0.87hk [7]. More typical thermal 
sources at ~ 2 jiK require velocity selection at the cost 
of atom number to achieve the required width. One can 
further reduce the momentum width of a condensate by 
manipulating trap parameters [13] or tuning the atomic 
interactions using a Feshbach resonance. BECs also have 
a transverse spatial width over an order of magnitude 
lower than a typical thermal source. Coupled with a low 
transverse momentum width, this leads to relaxed con- 
straints on spatial wavefront aberrations and intensity 
gradients in the LMT lasers. 

Our production of BEC is described in [19]. Briefly, 
we produce pure 87 Rb \F = 1, mp — — 1) condensates 
with up to 2 x 10 6 atoms in a crossed optical-dipole 
trap with cj x ,y,z = 27r x (50, 57, 28) Hz by evaporation 
in a magnetic and an optical trap. We then switch off 
the trap suddenly, allowing the cloud to expand and fall 
for up to 35 ms before probing the atoms using standard 
absorption imaging. We can transfer the atoms to the 
magnetically-insensitive \rriF = 0) state using a Landau- 
Zener rf sweep after the BEC is formed; however, we 
presently find this step unnecessary as we observe no ef- 
fect on fringe visibility, signal-to-noise ratio or our mea- 
surement of g. To produce a thermal cloud for compari- 
son, we terminate the magnetic trap evaporation earlier, 
loading fewer atoms into the optical trap. The final evap- 
oration then results in a phase space density just below 
that required for condensation at 100 nK. After releasing 
the cloud we wait 12 ms before initiating our interferom- 
eter cycle. 

Building a gravimeter using a BEC source allows spa- 
tially resolved imaging of the different momentum states, 
eliminating the need for state labelling via Raman transi- 
tions. This cancels systematic shifts [16] and greatly sim- 
plifies our laser system. We derive the two phase-locked 
optical frequencies required to drive Bragg transitions 
from an amplified external-cavity diode laser giving 1.3 W 
of light red detuned 90 GHz from the \5S 1 / 2 ,F = l) -» 
|5P 3 / 2 ,i^ / = 2) transition. The light is split through two 
80 MHz acousto-optic modulators (AOMs) driven by a 
direct digital synthesizer (DDS) referenced to a rubid- 
ium frequency standard. This allows us to produce the 
arbitrary pulse shapes and coherent frequency sweeps re- 
quired for the experiment. The first order of each AOM is 
combined on a polarizing beamsplitter and coupled into 



a polarization-maintaining fiber, resulting in 150 mW in 
a collimated 3 mm beam directed vertically through our 
science cell. The beam passes through a A/4 wave- 
plate before being retro-reflected by a mirror mounted 
on a multi-layer passive vibration isolation system. This 
scheme also generates a second pair of Bragg frequen- 
cies, but as we allow 12 ms of free fall to reach ballistic 
expansion before initiating the interferometer sequence, 
this pair is Doppler shifted well off resonance. We mea- 
sure the relative phase noise of the Bragg laser system to 
be negligible compared to that introduced by our retro- 
mirror. 

The Bragg resonance condition is given by S n = 4ncj r , 
where S n is the frequency difference of the Bragg beams 
for an nth order transition, and u r = hk 2 /2m is the 
single-photon recoil frequency. To operate in the Bragg 
regime and address the entire cloud, we use Gaussian- 
shaped pulse envelopes [20] and choose our pulse length 
r to satisfy ^ • k < r _1 < cj r , ensuring minimal loss 
to adjacent momentum states for a given n. For n = 1, 
we are able to maximise our 7r-pulse efficiency to 95% in 
this way. However, we find that a 300 /is velocity selec- 
tion pulse is required in order to achieve 93% efficiency 
for n = 3, an observation that supports the argument 
given earlier that the requirement on Ap becomes more 
stringent with n. Due to the size of our science cell, we 
are limited to an interrogation time of T = 5 ms. After 
several milliseconds of further separation following the 
final 7r/2 beamsplitter, an absorption image is taken to 
determine the number of atoms in each momentum state. 

The freely-falling atoms experience a time dependent 
Doppler shift Sd(t) = 27raot where ctQ = ^k • g is a fre- 
quency chirp. This modifies the Bragg resonance condi- 
tion in the laboratory frame to S n (t) = Anuo r J r2k • gt. We 
compensate this by sweeping 5 at a rate a c± 25.1 MHz/s, 
determined by local g near our lab in Canberra, Aus- 
tralia [21 . The interferometer phase then becomes = 
n(2k • gT 2 — 2iiaT 2 ). By scanning the sweep rate a, we 
record interference fringes with a period of 1/nT 2 . 



III. GRAVIMETRY WITH BEC 

Fig. [I] shows gravimeter fringes for n = 1, T = 3 ms. 
We observe a high visibility of (83 ±6)%. Increasing 
T generally reduces the visibility, and more rapidly for 
larger n. We speculate that wavefront aberrations in the 
Bragg laser beams contribute to this. Aberrations cause 
different atomic trajectories to experience different phase 
shifts, as has been discussed in [22j[23]. These different 
phase shifts are averaged through detection, causing a re- 
duction in fringe visibility. This effect will be exacerbated 
for larger n and T, as each atomic trajectory samples 
more of the transverse phase profile of the beam. We ex- 
pect significant aberrations in the current apparatus due 
to the close proximity of the Bragg beams to our mag- 
netic trapping coils, and suspect that with closer to ideal 
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FIG. 1: Interference fringes from a BEC-based gravimeter 
with n = 1, T = 3 ms. We observe a visibility of (83 ± 6)%. 
The solid line is a least-squares sinusoidal fit to the data (as 
in all figures). 



optical wavefronts, visibility would scale more weakly 
with n and T. Nevertheless, we are able to improve our 
sensitivity by scaling the interferometer space-time area 
using third-order Bragg LMT beamsplitters. We achieve 
a mid- fringe precision of A<£/<I> = 5 x 10 _4 Hz _1 ^ 2 in this 
way. The corresponding fringes are given in Fig. |2|a). 
These data represent 16 minutes of acquisition time. We 
can determine gravity from ao to be g = 9.7859(2) ms -2 . 
In [21], g is measured to be 9.795499189(29) ms" 2 , ap- 
proximately 11 km from our lab and 150 m higher in ele- 
vation. Our value disagrees at the 10 -2 level. This low 
accuracy is almost certainly due to the alignment of our 
Bragg beam along g, as when calculating g from ao, we 
have assumed k ■ g = kg. We estimate an alignment un- 
certainty of 3° in the current apparatus, which leads to 
a systematic error in g of up to 0.026 ms -2 . 

It is worth noting that state-of-the-art (SO A) BEC ma- 
chines can produce 2.4 x 10 6 condensed atoms/s [24 . The 
best published state- and velocity-selected thermal cloud 
used in a gravimeter has a factor of 25 higher flux at a 
momentum width of 0.87hk [7 j. It may be the case, how- 
ever, that when utilising increasingly large momentum 
transfer, the high spectral density of condensed sources 
will result in a higher usable flux due to the strict re- 
quirements on Ap, as discussed earlier. We are currently 
working on quantifying this effect in detail, and this will 
be the central topic of an upcoming paper. 



IV. THE EFFECT OF ATOMIC 
INTERACTIONS 

Fig. [2^b) shows the atom number for each correspond- 
ing point in the fringes in Fig. [5|a). The variation in 
number at each point is intentionally imposed. Despite a 
variation of 300% in density, we observe no detrimental 
effect of dephasing on signal-to-noise, or our measured 
value of g at our limit of precision. After 12 ms, the 
momentum width of the cloud is within 1% of its asymp- 
totic value, as shown in Fig. [5|c). We can estimate the 
interaction-induced phase uncertainty for our or a simi- 
lar device using the following simple model. For a 50/50 
beamsplitter on N atoms, the variance in the number 
difference for the two output modes is ~ y/N, giving 
a variance in the density of each mode. As a result of 
the mean- field energy shift, there will be an uncertainty 
in the relative phase evolution due to the energy uncer- 
tainty: 
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where U = Airh 2 a/m is the interaction parameter from 
the Gross-Pitaevskii equation, which is directly propor- 
tional to the mean-field energy, with a the scattering 
length. The initial peak atom density is n(0), and V(t) 
represents the cloud volume during expansion from a 
harmonic trap, which we can calculate from the Gross- 
Pitaevskii equation. Thus the uncertainty in the relative 
phase evolution is: 
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where \i ~ n(0)U is the chemical potential. This dephas- 
ing rate is then integrated through a given interferometer 
sequence to determine the phase uncertainty due to in- 
teractions, A<£ m f. It is worth noting that as this estimate 
assumes a uniform initial density equal to the peak den- 
sity in trap, it overestimates the effect. 

For our current experimental parameters we find that 
interaction-induced dephasing would limit precision to 
10 -7 per shot, well below our current sensitivity. Fig. 
|2|d) projects this estimate towards SOA device parame- 
ters, plotting the dephasing-limited sensitivity as a func- 
tion of T, for 10 6 condensed atoms/s and our trap pa- 
rameters. The shaded region represents the dephasing- 
limited sensitivity for expansion times ranging from 
^exp = 12 — 40 ms from the right to left boundary. The 
solid curve is the shot-noise limited sensitivity, and the 
dashed line the value of the current SOA sensitivity for 
an atomic gravimeter with T = 0.4 s and 6 x 10 T atoms 
[7]. We find that with an appropriate choice of trap 
parameters and expansion time, interact ion- induced de- 
phasing can be made negligible compared with the shot 
noise limit. More importantly, as T increases dephasing 
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quickly becomes negligible compared with the current 
SOA precision. Very recent work has comprehensively 
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FIG. 2: (Color online) Density related dephasing effects, (a) 
Fringes for our highest precision configuration of n — 3, 
T — 4ms. (b) Atom number in each run of (a), (c) Ver- 
tical momentum width of the BEC as a function of expansion 
time. Data points are measured using Bragg spectroscopy as 
shown in the inset. The solid lines are from a numerical sim- 
ulation of the Gross-Pitaevskii equation for our experimen- 
tal setup with no free parameters, (d) Estimated dephasing 
limited sensitivity as a function of interrogation time. The 
shaded region represents values for a range of expansion times 
(texp = 12 — 40ms from right to left). The solid line is the 
shot-noise limit for 10 6 atoms. The dashed line is indicative 
of the current SOA for an atomic gravimeter (with n = 1, 
T — 0.4 s), highlighting that for a significant range of pa- 
rameters, dephasing will not limit sensitivity compared with 
current SOA. 



investigated the effects of atom interactions in free-space 
BEC interferometers, and also confirm that interaction 
can be made negligible [25] . 

V. COMPARING A THERMAL AND 
CONDENSED SOURCE 

In Fig. [3j we give a comparison of fringes using a BEC 
and a 100 nK (T/T c > 1) thermal state as the source for 
our gravimeter, where we make every effort to ensure that 
the system is otherwise identical. In particular, we use an 
identical velocity selection pulse for each sequence. The 
fringe data sets were taken consecutively. The condensed 
source shows an improved result compared to the ther- 
mal state, with the visibility increasing from (58 ±4)% to 
(85 ± 11)%. A 500 nK thermal state gives an even lower 
visibility of (50 ±5)%. As discussed earlier, we speculate 
that wavefront aberrations contribute to the observed dif- 
ference, as at 100 nK the thermal cloud has a factor of 3 
larger transverse momentum width than the condensate. 
Thus it will have wider range of atomic trajectories, and 
therefore a wider range of phase shifts across the cloud 
[23] . This observation suggests that even a factor of 3 
smaller transverse momentum width can improve fringe 
visibility if aberrations limit the interrogation time. 
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FIG. 3: Comparison of thermal and Bose-condensed atomic 
sources. We use a n — 1, T — 3 ms gravimeter cycle. Both 
have an identical vertical momentum width, but differ in their 
transverse momentum width. A significant improvement in 
visibility from to 58% to 85% is seen for the BEC, with all 
other experimental parameters kept constant. 



VI. BLOCH OSCILLATION BASED LMT 

We also achieve a comparatively high fringe visibility 
of (24 ± 4)% with 6Hk LMT beamsplitters using Bloch 
oscillations (Fig. |4|l). After an initial 4hk Bragg 
beamsplitter, we adiabatically load a lattice of depth 
~ 10E r with q = (stationary in the atom frame) 
in 100 /is, where q is the quasimomentum and E r the 
single-photon recoil energy. We then chirp 5 over 200 /is, 
sweeping q through one Brillouin zone. The momentum 
of one arm is thus increased by 2hk in the lab frame, 
as one arm remains in the lowest band whilst the other 
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FIG. 4: (a) Fringes from a LMT gravimeter using Bloch beam- 
splitters. A visibility of 24% is observed for T = 2.5 ms, 
and effective order n = 2.42, calculated from the space-time 
area. This is in agreement with the fitted fringe period of 
(70 =b 5) kHz. (b) The intensity of the pulse sequence used 
for this interferometer, and the resulting space-time diagram. 
Using only the Gaussian pulses results in a standard Mach- 
Zehnder interferometer, (c) Absorption images showing the 
two arms of the interferometer after each pulse. The scale bar 
in the (c) represents 300 /im. 



is increased beyond !0E r , or our acceleration time is 
increased to impart larger momenta, the interferometer 
output converges to P = 0.5. We are currently inves- 
tigating this effect further, and do not believe it to be 
due to randomized light-shift-induced phases as these 
would tend to reduce visibility without loss of contrast. 
We have been able to apply a 30hk beamsplitter to our 
BEC with an efficiency > 95%, limited so far by the size 
of our science cell. 

VII. CONCLUSIONS 



The question of whether to use coherent or thermal 
atomic sources for precision inertial sensing remains an 
important one requiring further investigation. We have 
presented results from the first comparison of these in 
a Mach-Zehnder gravimeter. We observe interference 
fringes with a high visibility for the condensed source, 
and are able to increase our sensitivity to gravity by im- 
parting larger momentum to the atoms in the beam split- 
ting process. The thermal source produces fringes with 
significantly lower visibility than the condensed source. 
We believe this is a result of its larger transverse mo- 
mentum width, which causes higher sensitivity to wave- 
front aberrations in the beamsplitter lasers. We have 
also presented a simple model which demonstrates that 
interaction-induced dephasing is negligible in an interfer- 
ometric measurement with freely-falling, coherent atomic 
samples. It may be the case that exploiting the sub-recoil 
momentum distribution of Bose-Einstein condensates to 
realise high visibility fringes with very large momentum 
transfer beamsplitters will lead to sensitivity beyond cur- 
rent state-of-the-art in precision inertial sensors. 



undergoes inter-band Landau- Zener transitions [TT] . 
This process is reversed to decelerate this arm before 
the 7r-pulse, after which the other arm is subjected to 
the same procedure (Fig. [4)3). We use T — 2.5 ms, and 
our pulse sequence gives the interferometer a space-time 
area with an effective n = 2.42, in agreement with the 
fitted fringe period. In contrast to the work in [lOj [TT] , 
we have not used symmetric acceleration of each arm to 
balance differential light shifts, although our space-time 
area is smaller. We find that if either the lattice depth 
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